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ABSTRACT 



We have reviewed the chemistry and cooling behaviour of low-density (n < 10 4 cm" 3 ) primor- 
dial gas and devised a cooling model wich involves 19 collisional and 9 radiative processes and is 
applicable for temperatures in the range (1K< T < 10 8 K). We derived new fits of rate coefficients 
for the photo-attachment of neutral hydrogen, the formation of molecular hydrogen via H~ , charge 
exchange between H2 and H + , electron detachment of H~ by neutral hydrogen, dissociative recom- 
bination of with slow electrons, photodissociation of H^, and photodissociation of H2. Further 
it was found that the molecular hydrogen produced through the gas-phase processes, + H — > 
H2 + H + , and H + H — > H2 + e~ , is likely to be converted into its para configuration on a faster 
time scale than the formation time scale. We have tested the model extensively and shown it to 
agree well with former studies. We further studied the chemical kinetics in great detail and devised 
a minimal model which is substantially simpler than the full reaction network but predicts correct 
abundances. This minimal model shows convincingly that 12 collisional processes are sufficient to 
model the H, He, H + , H , He + , He ++ , and H2 abundances in low density primordial gas for 
applications with no radiation fields. 



1. Introduction 



Since Szalay (1967) realized the importance of the gas phase reactions for H2 formation in primordial gas 

H + e" -» H + hu, (1) 

H + H -> H 2 + e" (2) 

and 

H+ + H -> H+ + hi/, (3) 

H+ + H -> H 2 + H+, (4) 

and Peebles and Dicke (1968) formulated their theory for globular cluster formation based on H2 cooling of 
collapsing primordial density fluctuations, many models have been put forward employing this effect to explain 
the origin of structures on a vast range of mass scales. E.g., for the cosmic string model of structure formation, 
molecular hydrogen is believed to trigger the fragmentation of late cosmic string wakes. This is because, for the 
weak accretion shocks encountered during the formation of late cosmic string wakes, thermal instabilities due to 
hydrogen line, or bremsstrahlung cooling are not accessible (Rees 1986). Kashlinsky and Rees (1983) discussed 
fragmentation of primordial gas that cools via H2 . Couchman and Rees (1986) envisaged star cluster sized 
objects collapsing via H2 cooling at high redshift which then reionize the intergalactic medium. More recently, 
a model for the origin of halo globular clusters and spheroid stars based on a cooling instability triggered by 
H2 formation has been put forward by Vietri and Pesce (1995). Furthermore, in studies of primordial star 
formation, molecular hydrogen cooling plays a central role (see Stahler 1986 for review). All these applications 
are of fundamental importance for our understanding of the origin of structure in the universe and its subsequent 
evolution. 

All of the above problems are inherently multi-dimensional and a method that computes the hydrodynamics 
along with the chemistry is desirable. In order to form H2 efficiently through the above gas phase reactions, 
electrons and protons have to be abundant at relatively low temperatures (T < 10 4 K) allowing the formed 
molecules to survive. This situation arises naturally in shock heated gas, where the gas recombines slower than it 
can cool, and an enhanced ionized fraction (over the equilibrium value) is reached, despite the low temperatures. 



This has been convincingly shown by many previous investigations (Hollenbach and McKee 1979, MacLow and 
Shull 1986, Shapiro and Kang 1987, Anninos and Norman 1996). Furthermore, at recombination the universe 
expands so fast that recombination will not be complete and the gas is left with a residual abundance of free 
electrons (see Peebles 1993). This residual ionized fraction can be enough to form a substantial amount of 
H2 in structures that only have very weak virialization shocks (see Couchman and Rees 1986; Tegmark et al. 
1996). However, the chemistry has to be solved self-consistently with the structure formation equations, which 
is computationally due to the stiff nature of the reaction network. This has forced former studies to constrain 
themselves to zero dimensions, using e.g. steady state shock approximations. Only recently Haiman et al. 
presented a study of collapsing small scale structure that incoorporated the time dependent chemistry in a one 
dimensional hydrodynamics code. In our methodology paper (Anninos, Zhang, Abel, & Norman 1996, hereafter 
AZAN96) we discuss a numerical method that unifies the time dependent chemistry with multi-dimensional 
cosmological hydrodynamics. With the chemistry model which accurately predicts abundances, and the cooling 
behaviour of primordial gas for temperatures from one to ~ 10 8 Kelvin and all densities below ~ 10 4 cm -3 , this 
introduces a powerful tool in investigating many aspect of structure formation by means of three dimensional 
hydro dynamical simulations. 

The paper is organized in the following way: In the next section we review some general properties of 
primordial gas such as nucleosynthesis constraints and the coronal limit. In section 2.2. we discuss general 
arguments on how to select the dominant reactions. We then give an elaborate presentation of the collisional 
and radiative processes which we find to be important and the rate coefficients we have adopted. Section 5. 
gives an overview of the processes we find to be negligible for the range of applications we are interested in. 
Then we review all cooling mechanisms needed for our model and discuss the molecular cooling and heating 
functions in the low-density limit. Finally we present an extensive discussion of the performance of our model 
for an application to a strong shock wave arising from the collapse of a single pancake or cosmological sheet. 
The complete chemical model is summarized in appendix A.. The discussion why we find all H2 to be in its 
para configuration is given in appendix B.. 



What are the major species determining the physics in primordial gas? Here we briefly summarize what we 
include in our cooling model and defer the detailed discussion to the following sections. 

We know from nucleosynthesis and observational constraints that 7 Li/H ~ 1CT 10 , D/H~ 10 -5 , 3 He/H ~ 10~ 5 , 
and 0.236 < 4 He/H < 0.254. Where H, D, He, Li denote the mass abundances of hydrogen, deuterium, helium 
and lithium respectively. Hence D, 3 He, and Li have, compared to neutral hydrogen, very low abundances. For 
the line cooling of D, 3 He, and Li to contribute significantly to the overall cooling of the gas, their excitation 
rates would have to be many orders of magnitude higher than the corresponding processes of neutral hydrogen. 
However, this is not the case and those species, at least in their atomic and ionic forms, will have negligible 
influence on the hydrodynamics of the gas. The work by Lepp and Shull (1984) indicates that molecular 
hydrogen was by far the most abundant molecule in the universe at redshifts between recombination and before 
the first stars formed. It forms by radiative association of the negative hydrogen ion with neutral hydrogen 
and charge exchange between H and Hj. The minimal temperature at which rot./vib. of molecules can get 
excited defines a minimal temperature, and hence a minimal Jeans mass. In the case of molecular hydrogen, 
this introduces a minimum temperature of about 10 2 K (see also Mac Low and Shull 1986, Shapiro and Kang 
1987) which translates to a Jeans mass of 



2. General Properties of Primordial Gas 
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which tells us that in an application as e.g. proto-galaxies, where we expect number densities of the order 
100cm -3 we expect a Jeans Mass of about 10 4 M Q . This is a rough estimate of the minimum mass scale one has 
to resolve in numerical simulations. If one includes also HD and LiH which will get excited only for temperatures 
above 112K and 21K, respectively, one has be able to resolve masses nearly a magnitude smaller. A study by D. 
Puy et al. (1993) was able to give fairly good estimates on the primordial HD abundance after recombination. 
They found that HD was the second most abundant molecule with uhd ~ 10~ 3 ' 5 n# 2 . We know, therefore, 
that molecular hydrogen will always be the dominant molecular coolant for temperatures above ~ 100K. Hence 
choosing to not include other molecules than H2 can be understood that we restrict our attention to masses 
greater than ~ lO 4 M0(n/lOOcm -3 ). Recently the Lithium Chemistry has been rigorously discussed by Stancil 
et al. 1996, where they found that much less primordial LiH is formed than previously expected and suggested 
that a significant fraction of LiH might only be formed by three body reactions. Since our model focuses at 
number densities below ~ 10 4 cm -3 and three-body reactions typically become important at greater densities 
we conclude that in the density regime we are interested in LiH will not be important. The papers mentioned 
above did not include H-j~ which, due to its low abundance is not believed to be a significant coolant. From the 
calculations by Lenzuni, Chernoff, and Salpeter (1991) we know that in fact the H3" abundance can be of the 
same order as the H^~ abundance, but their results also show that it has a negligible effect on the molecular 
hydrogen abundance. We also did not include He molecular ions as well as hydrogen ion clusters (H+ with 
n > 4) because their expected abundances are even smaller than the ones of the molecules mentioned above. 



2.1. "Molecular" and "Atomic" Coronal Limit 

Atoms and molecules have very complex intrinsic properties, and their chemical behavior varies somtimes 
drastically with the specific quantum-mechanical state they occupy. For atoms and ions at moderate or low 
densities like in the solar corona (electron number density n e ~ 10 s — 10 9 cm~ 3 ), the following features of 
thermodynamical equilibrium do not hold (Sobelman et al. 1979): 

• Boltzmann distribution of atoms over excited states. 

• Saha distribution of atoms over degrees of ionization. 

• Principle of detailed balance. 

The velocity distribution of the free electrons is, however, as a rule nearly always Maxwellian. In this 
low-density limit we know that the level distributions are given by, 

/ bp 

N- A k 

where Nf denotes the number density of the species i in its level k, vo~jk the rate coefficient for collisional 
excitation from level j up to level k, the total probability for spontaneous transition from all higher levels 
down to k. This approximation is applicable, if, 



n e < 



{vo-jk) 



One important assumption here is that collisional excitations outweigh radiative excitations which is always 
true as long as there are only moderate external radiation fluxes. For the s-p excitation in Helium one finds 



4 



that the coronal limit holds up to electron densities of ~ 10 17 cm 3 . ( Here we used {vajk) ~ 6 x 10 9 cm 3 s 1 
at 400K from Janev et al. 1987, (2.3.1), and A k = 1.8 x 10 9 s" 1 (Sobelman et al. 1979, p297)). Thus the 
coronal limit is valid for will be applicable in all our intended applications concerning structure formation. 
In proto-galactic clouds for example, we expect total number densities of about 100 cm -3 . For us the most 
important point is that we find nearly every atom in its ground state. Thus it is not necessary to treat multilevel 
atoms in our case. 

Comparing the time scales for collisional excitation, collisional de-excitation (Lepp and Shull, 1983) and also 
the transition probabilities (Allison and Dalgarno, 1970) of molecular hydrogen, one finds that for low densities 
{rin < 10 4 cm -3 the population of excited states of H2 is by many orders of magnitude smaller, than the ground 
state. This is a fact which can also be seen in Lepp and Shull 1983, as well as in the work of Dove et al. 1987 
where they studied dissociation of molecular para-hydrogen by collisions with helium. They found that at low 
densities the dissociation rate coefficient approaches a constant value, which corresponds to direct collisional 
dissociation out of the (v = 0, J = 0) level only. A strong radiation field in the UV is capable of populating 
the exciting states and hence changing the chemical behaviour of H2 . However, a detailed study (Shull 1978) 
shows that only for fluxes greater than ~5x 10 -15 erg cm~ 2 s _1 Hz~ x sr _1 UV Pumping becomes effective. 
We conclude that in our intended calculations we can assume all our considered species to be in their ground 
state. 

2.2. Selection of Reactions and Rate Coefficients 

Although primordial gas is a simple mixture of hydrogen and helium, the ongoing physical reactions in it are 
immense. E.g. in Janev et al. 1987 one finds more than 70 reactions only involving the ground states of our 
species. It is obvious that in order to construct a computationally feasible model that some selections have to 
be made. The quest is to make them as good as possible i.e. without neglecting any important physics. To 
neglect a collisional processes one, obviously, has to make sure that 

1. for all reactants and products the rate will for any temperature and density never dominate the right 
hand side of the rate equations. 

2. the reaction enthalpy is always negligible energy contribution or loss to the internal energy of the gas. 

We used these criteria to construct our reaction network out of hundreds of reactions found in the literature as 
well as in databases. 

3. Collisional Processes 

In appendix A., we present all included processes, their rate coefficients and the according reference. 
Discussing the reliability of the used rate coefficients and their sources will be the main aim of this sec- 
tion. Further information on atomic rates can be found at Dima Verners 1996 Atomic Data Page 
(http : / / www . pa . uky . edu/ verner/ atom . html) 

3.1. Ionizing processes 

For the cosmological problems we are interested in, we have seen above that spontaneous decay is a lot faster 
than collisional excitation (coronal limit). This already pointed out the non-LTE character of our problem. 
Metaphorically speaking, the radiation field is in our case much smaller than in the LTE case, since there are 
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not enough collisional excitations to build up a Planckian spectrum. Therefore, if we compare in the follow- 
ing collisional and radiative processes, we can use arguments derived for strict LTE for qualitative statements 
on non-LTE conditions by assuming that the radiation temperature is much less than the matter tempera- 
ture. This is always true in the matter dominated phase of the universe, when radiation fields other than the 
CBR arc negligible. We adopt here Mihalas' (1978, pl23) estimate of the ratio between the total number of 
photoionizations Rik and the total number of collisional ionizations Cjfe, 



R ik 4(2vr 3 )^ 3 / Wk B T R \ 

~Ff — ~ 1 r ex P 

Ufc 3ni2e 2 h 2 c 3 \n e (k B T e )2 J 
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k B T e k B T R 



e 



where E{ is the ionization energy out of the level i, k B is the Boltzmann constant, m is the electron mass, T t 
the electron temperature, e the electron charge, c the speed of light, and W is defined by J v = WB U (T R ), where 
J u is the specific intensity at the frequency v, and B u (Tji) denotes the Planck spectrum of a black body with 
temperature, T R . 

This gives for ionization of hydrogen out of the ground state (Ei = 13.6eV) with the number density of free 
electrons in cm -3 , 



Rik iaW ( Tl\ 2 15782o( ±- 

~ 1.55 x 10 14 — -# e V Te 



Cik n e \ T ( 



e 



In thermal equilibrium with T e = T R ~ 10 4 and n e ~ 1cm 3 this ratio is of the order 10 15 , therefore radiative 
ionization dominates by far. For T e ~ 10 x T R ~ 10 3 , however, the ratio 5^ is of the order 10 -47 and we 
see that for all of our intended applications collisional ionization is the dominating process. It is important, 
however, to include photoionization to study the influence of external radiation fields on the chemistry and the 
cooling to heating balance of the structure forming gas, as, e.g., Lyman Alpha clouds in the vicinity of quasars 
which are observationally accessible. 



(1) Collisional Ionization of Hydrogen We use the fit of Janev et al. (1987, 2.1.5). For temperatures 
below 1000K~ 0.086eV we take it to be zero since it is smaller than 10 _60 cm 3 s _1 . Since in our simulations 
we do not expect to encounter temperatures above 10 5 eV no relativistic corrections were made. The ionization 
threshold is 13.6eV. 



(2) Collisional Ionization of Helium We use the fit given by Janev et al. (1987, 2.3.9) For temperatures 
below 4000K we take it to be zero since it is smaller than 10 -38 cm 3 s -1 . The threshold lies at 24.6eV. 



(3) Collisional Ionization of He + The rate is given in the Aladdin Database (1989) of the IEADS (In- 
ternational Atomic Energy Agency, Data Section). The energy threshold is 54.4eV. Again we only have to 
consider ionization out of the ground state and no relativistic effects. 



3.2. Radiative Recombination 



Radiative recombination is the inverse reaction of photoionization. With a similar argument as the one we 
used to discuss the relative importance of collisional to radiative ionization, one finds (Mihalas 1978) that 
radiative recombination always outweighs the collisional one. This is especially valid in our low-density limit 
since the collisional recombination which is, as all three-body processes, typically negligible for densities smaller 
than about 10 s cm -3 . 
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(4) Radiative Recombination to Hydrogen In the coronal limit we assume that recombination can 
happen into any quantum state which, if it is an excited state, will spontaneously decay to the ground state. 
Therefore, the total rate coefficient is the sum of the rate coefficients a n for all n = I..00. Ferland et al. (1992) 
computed hydrogenic rate coefficients which are in principle exact. We have fitted their data (the sum of all 
rate coefficients for n = 1..1000) to a form similar to what Janev et al. (1987) used in their compendium. We 
made sure that the fit is accurate for temperatures from IK to 10 9 K. 

(5, 6) Recombination to Helium He + is the only species in our model, which is subject to di-electronic 
recombination which dominates at high temperatures (T > 6 x 10 4 K). Since radiative and di-electronic re- 
combination rates are independent of density their sum gives the overall recombination rate. For the radiative 
recombination process we employ the rate coefficient given by Cen (1992) and for the di-electronic recombination 
the one given by Aldrovandi & Pequignot (1973). 

(7) Photo-attachment of Hydrogen The rate coefficient from Hutchins (1976) is stated to be accurate to 
10% in the temperature range 100K(= 0.0086eV) < T < 2500K(= 0.254eV). The cross section for the inverse 
reaction has been calculated by Wishart (1979) to within 1% around the threshold (2 x 10 14 Hz< v < 2 x 10 15 Hz). 
From that data alone one can compute the rate for photo-attachment from 2000K to 10000K, using the principle 
of detailed balance or simply the Saha equation. To cover a greater temperature range, however, we use the 
cross section given in de Jong (1972) for all frequencies outside of the interval given by Wishart (1979). With 
that we are able to compute and fit this rate in the temperature range from IK to 10 8 K, with an accuracy to 
within a few percent for IK < T < 100K, about one percent for 100K< T < 10 4 K, and better than 10% for 
T > 10 4 K. At LTE the Saha equation has to be valid and the rate at non-LTE will be naturally the same as in 
LTE since it is an atomic property. The Saha-Boltzmann ionization equation reads, 

kdet _ Ny+i „ _ Ur+i 2(27rmfcT) 3 / 2 1^ 

k att N r 6 U r h* ^ kT*' 

where N e denotes the free electron number density, N r+ i the number density of neutral hydrogen atoms, N r 
the number denisty of H~ ions, kdet the rate coefficient for photo-detachment, k a tt the rate coefficient for photo- 
attachment, u r , u r+ \ the respective partition functions, I r the threshold for photo-detachment of H (0.755eV), 
m the electron mass, and k the Boltzmann constant. The single bound state of H is l S and the ground state 
of H is 2 S. Hence u r +i = 2 and u r = 1. The rate coefficient, kdet, for photo-attachment in LTE is derived by 
integrating the cross section over a Planckian black body spectrum. The fit, which is given in the appendix, 
equals the one given by Hutchins (1976) to within a few percent in the range where the latter is applicable. 

(8) Formation of molecular Hydrogen via H The rate coefficient has been taken from Janev et al. 
(1987, 7.3.2.b) and is based on theoretical cross sections of Browne and Dalgarno (1969) and normalized to the 
experimental results of Hummer et al. (1960) in the region of 500eV to 4 x 10 4 eV. We averaged their energy 
dependent rate coefficient over a Maxwellian velocity distribution of the incident particle. The used data is 
only reliable for temperatures above O.leV. However, we, know that the rate coefficient is constant for small 
temperatures and we also see that the energy dependent rate coefficient is already the same for O.leV and 
leV of the incident particle for a temperature of O.leV. Therefore, we take the rate coefficient at O.OleV to 
be the same as at O.leV, fit it, and set it to be constant for lower values. This leads to a value of the rate 
coefficient for low temperatures about 10% higher than the constant values given by Bieniek (1980) and de Jong 
(1972). Launay, Le Dourneuf, and Zeippen (1991) computed the ratecoefficient for this reaction depdening on 
the ro-vibrational state of the produced H2 molecule. Using the potential derived by Senekowitsch et al. they 
found a thermal rate at 300 K of 1.49 x 10~ 9 which lies 4% above the rate we derived from the data of Janev 
et al. (1987). 
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(9) Formation of The rate coefficient for this radiative association reaction has been calculated by 
Ramaker and Peek (1976) and has then been fitted by Shapiro and Kang (1987) as well as by Rawlings et al. 
(1993). We prefer the fit given by Shapiro and Kang, because it covers a greater temperature interval. 

(10) H2 Formation via The rate constant for this charge exchange reaction between neutral hydrogen 
and H^ was measured by Karpas et al. (1979) to an accuracy better than 20%. This rate is constant at low 
energies because the distribution of final states in phase space is determined by the energy released in the 
reaction, which is almost independent of the incident energy (Peebles 1993). By assuming it to be constant at 
high temperatures as well, we potentially introduce an error in the H2 abundance, because the by far dominating 
destroying processes will be balanced by a wrong formation rate. However, this error will not be significant 
because the formation of hydrogen molecules is dominated by the H formation path. Figure 1 summarizes 
the rate coefficients for formation of H and as well as the ones for H2 . 

(11) Charge Exchange between Molecular Hydrogen and H + This is the inverse reaction of the above. 
The rate coefficient was derived by Donahue and Shull using detailed balance and the Karpas et al. (1979) data 
to 

10 / 2.65eV\ 
k u =6.4 x l(T 10 exp ^ ^—j . 

That rate differs drastically from the one we have computed from data given in Janev et al. (1987) which is 
accurate also for high temperatures. The data of Janev et al. has to be considered by far more accurate than 
the simple detailed balance argument of Donahue and Shull (1991). However, comparing results of our model 
with the one of Shapiro and Kang (1987) we find that drastic differences in this rate do not change the final 
abundances significantly. 

(12) Dissociation of molecular Hydrogen by Electrons We employ the rate coefficient that Donahue 
and Shull (1991) derived from the cross section given by Corrigan (1965). 

(13) Dissociation of Molecular Hydrogen by neutral Hydrogen Dove and Mandy used fully 3D quasi- 
classical trajectory calculations, using the Siegbahn-Liu-Truhlar-Horowitz potential energy surface, to calculate 
the rate coefficient in the low density limit. They found significant differences to the former "standard" values 
given by Lepp and Shull (1983). We use the results for the total rate coefficient of Dove and Mandy (1986), 
which assumes that all the quasi-bound states will dissociate. Although their ratecoefficient is not accurate at 
temperatures well below 500 K this does not lead to a significant error in the H2 abundance as can clearly be 
seen from the discussion in section 7.1. 

(14) Collisional Detachment of H by Electrons We take the rate coefficient from Janev et al. 1987 
(7.1.1), which is based on experimental results. The net energy loss is as in ionization reactions the threshold 
energy which is 0.755eV. 

(15) Collisional Electron Detachment of H by neutral Hydrogen The formation of an auto- 
dissociating state of molecular hydrogen, (equation 13 in section 5.6.) seems to be far more effective in destroying 
H~ than the direct collisional detachment 

H + H -> H + H + e~. 

The rate coefficient for the auto-dissociating state below O.leV is very uncertain. We have extrapolated the 
data given by Janev et al. (1987) linearly for those temperatures. This overestimates the rate coefficient in 
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this range. However, it does not introduce any significant error in the abundances since for a collision of H 
with H in this temperature regime it is far more probable to form a stable H2 molecule than to destroy the H~. 
In fact, we expect this reaction to be only of importance in a very small range around leV, since for higher 
temperatures, the dissociation by free electrons is far more effective than the one by H. 

(16) Mutual Neutralization between H and H + The rate coefficient given by Dalgarno and Lepp 
(1987) is consistent with the measurements by Szucs et al. (1984) and Peart, Bennett and Dolder (1985). See 
the latter reference for a review on the experimental difficulties and uncertainties of ion-ion collisions. 

(17) Formation in H _ and H + Collision Shapiro and Kang (1987) derived the rate coefficient from 
the cross section which was measured by Poulaert et al. (1978) in the range from O.OOleV to 3eV and found to 
have a E~ ' 9 dependence with 1.5 x 10~ 14 cm 2 at 0.003eV. The fit given in Shapiro and Kangs (1987), however, is 
discontinuous at 10 4 K. We changed it within 4%, which surely lies within the uncertainties of the experimental 
data to make it continuous. The result is given in the appendix. 

(18) Dissociative Recombination of with slow Electrons The rate for this reaction was calculated 
using multichannel quantum defect theory (MQDT) by Nakashima et al. (1987) in the energy range from 0.02 
to about leV. Their result is much smaller at low temperatures, than previously assumed (see e.g. Shapiro and 
Kang 1987, Mitchell and Deveau 1983). More recently Schneider et al. (1994) recalculated this rate coefficient 
also using MQDT. Their tabulated rate coefficient for dissociative recombination of the ground state via 
the E+ state is, above 1000K, by a factor of two to four smaller than the rough fit given by Nakashima. Our fit 
of the Schneider et al. (1994) data neglects the fact that the rate is increasing for T < 50K. This is reasonable 
since one does not expect this reaction to be important for temperatures smaller than 100K due the small 
abundance of free electrons, as well as H^. 

(19) Neutralization between and H This mutual neutralization reaction takes place at low energies 
(< leV) where H^ and H~ coexist. Unfortunately, the only three measurements we are aware of, Aberth et 
al. in Moseley et al. (1975), Szucs et al. (1983), and Dolder & Peart (1985), give only data for energies above 
3eV. We use the rate coefficient given by Dalgarno and Lepp (1987) despite its uncertainty. This equals the 
coefficient given by Prasad and Huntress (1980) which was used by Shapiro and Kang (1987) to within 25%. 



In general we will only consider photoionization (dissociation) out of the ground state since we assume the 
abundance of excited states to be negligible. The rate of a particular photoionization (dissociation) reaction is 
given by: 



where k denotes H, He + , H~, H2 and H2 respectively, i(v) is the intensity of the radiation field, v t h the threshold 
energy for which photo ionization(dissociation) is possible, and a* the frequency dependent photoionization 
(dissociation) cross section of species k. 



4. 



Photoionization and Dissociation Processes 




(7) 
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(20 — 22) Photoionization of H, He and He + Hydrogenic cross sections have been studied in great detail. 
Since they are relatively simple to calculate the available data is very accurate. We use the typical expression 
given in Osterbrook (1989). The threshold frequencies are hv = 13.6eV, 24.6eV, and 54.4eV, for H, He, and 
He + , respectively. 

(23) Photo-detachment of the H Ion The best data available is given by Wishart (1979) with a stated 
accuracy within 1%. The fitting formula for the cross section given in Shapiro and Kang (1987) is accu- 
rate to ~ 10% at the high energy end. This is a normal characteristic of such fitting formulas. Although 
is adequate for present purposes, we recommend using the tabulated values of Wishart (1979), where they 
are avaialabel, in order to derive the rate of H~ destruction since using the fitting formula in the integral 
of equ. 7 leads to an overestimation of the according rate. H~ gets destroyed efficiently by the CBR to 
redshifts ~ 110. A fit for ^23 depending on radiation temperature can be found in Tegmark et al. 1996 
(http : / / www . mpa-garching . mpg . de/ max/ minmass . html) . 

(24) Photoionization of Molecular Hydrogen The high threshold of 15.42eV assures that the CBR will 
not be able to photo-ionize H2 in the post-recombination universe. The rate is taken from Shapiro and Kang 
(1987). 

(25 26)Photodissociation of H^ 

There is a conceptual difficulty in the treatment of this process since for strong radiation fields and high 
densities, excited states of Hg" will be populated. For them the photodissociation cross sections are much 
higher than for the ground state. However, the levels might not be in their thermal equilibirum distribution. 
To treat this exact, one would have to include all H2" levels explicitly, which for our applications would not be 
worth the effort since most of the H2 is formed via H _ . By using solely the cross section for the ground state 
we will overestimate the H^ abundance. Using the LTE rate one could derive a lower limit and estimate the 
error made, which we did not attempt in our simulations, yet. The most recent data for dissociation out of 
the ground state as well as LTE rates can be found in Stancil (1994). We have fitted the cross section for the 
ground state to better than 2%. This fit is more accurate, and also does not show the unfortunate divergent 
character than the one given in Shapiro and Kang (1987) which was drawn from Dunn (1968). The threshold 
energy is 2.65eV. For photons with energies greater than 30eV, H^ can get dissociated with two protons as 
products (proc. 26). For this we adopted the rate given in Shapiro and Kang (1987). 

(27) Photodissociation of H2 by predissociation Photodissociation of the ground state of molecular 
hydrogen (X 1 S+(u)) happens mostly through absorption in the Lyman- Werner Bands to the electronically 
and vibrationaly excited states, -B 1 S+(v') and C 1 n u (w / ) which then decay to the continuum of the ground 
state. This is called the two-step Solomon process (cf. Stecher and Williams, 1967). Allison and Dalgarno 
(1970) computed the band oscillator strengths and Dalgarno and Stephens (1970) derived the probability that 
those states decay into the continuum of the ground state. The dominant photodissociation paths are through 
absorption in the Lyman Band to the vibrational states 6 < v' < 20. This means that H2 dissociation happens 
mostly in a very narrow frequency range \2.2AeV < hv < 13.51eV. Assuming the incident radiation field to be 
practically constant around hv = 12.87eV, which corresponds to v' = 13, and neglecting self-shielding we can 
derive a rate constant for photodissociation of H2 through 
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where j denotes the radiation flux in s _1 ergs cm _2 Hz _1 at hv = 12.87eV, /q' the oscillator strength for the 
transition X 1 ^^ = 0) to p" Ly the probability that i? 1 E+(w / ) decays to the continuum, <fi v '{v) the 

line profile, and ire 2 /(mc) = 2.65 x 10~ 2 cm 2 the classical total cross section. 



(28) Direct Photodissociation of H2 Direct photodissociation of the ground state by absorbtion into the 
continua of the Lyman and Werner systems has been studied by Allison and Dalgarno 1969. We have fitted 
their cross sections within the stated accuracy of their data with, 



where 



^ = £+t + <n + (i - + < l ) (8) 



CT^g = 10" 18 X I 



dex(15.1289 - 1.05139 x hi/) cm 2 , 14.675 eV < hv < 16.820 eV 
dex(-31.41 + 1.8042 x 10~ 2 (hz/) 3 - 4.2339 x 10~ 5 (/iz/) 5 ) cm 2 , (9) 

16.820 eV < hv < 17.6 eV, 

a wo = 10- 18 x { dex(13.5311 - 0.91826181uv) cm 2 , 14.675 eV < hv < 17.7 eV, (10) 

Li i 8 { dex(12.0218406 - 0.819429hi/) cm 2 , 14.159 eV < hv < 15.302 eV, ( 

0-28 '' ' X j dex(16.04644 - 1.082438hi/) cm 2 , 15.302 eV < hv < 17.2 eV, ^ ^ 

cr^ 1 = 10~ 18 x { dex(12.87367 - 0.85088597hi/) cm 2 , 14.159 eV < hv < 17.2 eV, , (12) 

and y is the ortho-, to para-H2 ratio. Since unfortunately they only gave the data for photon excess energies 
up to 3 eV one might make quite significant errors in deriving the ratecoefficient where one integrates over this 
cross section. The error made depends on the shape of the radiation spectrum. 



5. Brief Summary of Neglected Processes 



In this section we will briefly summarize some processes which we find to be negligible. Although we can 
almost be sure that this smmary will be far from complete, we hope to stimulate with this a discussion 
which should lead us to a more detailed understanding of the primordial gas chemistry and enable us to find 
the minimum set of reactions, which describes primordial gas accurately enough for the desired cosmological 
applications. 



5.1. Photoionization Heating and Secondary Ionization 



For all simulations where we do want to consider strong external radiation fluxes we have to look at the fate 
of the photo electrons produced in photoionizations. To estimate the importance of this effect we first derive 
the numbers of high energy photo electrons produced in a typical quasar radiation field. 

Shapiro and Kang (1987) modeled a quasar radiation flux with: 

F v = 1.0 x 10 _21 e (^-} er 9 s cm' 2 s' 1 Hz' 1 , 

where vh is the Lyman edge frequency, a = 0.7, 0.7eV < hv < 12.4keV, and e is a parameter introduced to 
adjust the amplitude of the flux (or effectively the distance of the quasar). Although Sargent, Steidel, and 
Boksenberg (1989) argued that a > 1 might be a better fit, we will derive an upper limit of produced photo 
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electrons with the more drastic value of a = 0.7. We calculated the fraction n low _ln %9 _ between the electrons 

7, e / 7, e 

produced by photons in the range from hvn to 2hvn and the ones produced by photons in the range of 2hvjj 
to 12.4keV: 

low r2u H (T20(v) J . 

"7, e~ . . Jvh u hv aV 

high ~ f 24keV/h „ cr 2u (V) , / 
"7, e~ hv H r v~hU~ aU 

For a very hard spectrum with a = 0.7 we get n l ° w e - /n^ ig ^_ = 10 whereas for a more soft spectrum this number 
increases (e.g. 13 for a = 1.0). Because photo electrons have an energy of h(y — uh), the ones produced by 
photons in the range from hvu to 2hvjj will not be capable of collisionally ionizing any species other than 
H . Therefore, all they can do is to either go into heat, which means scatter with other electrons or ions and 
equilibrate to a Maxwellian distribution, or excite other atoms, ions or molecules. The photo electrons above 
2huff carry enough energy to also ionize (dissociate) other species. The very high energetic ones are capable of 
ionizing (dissociating) many atoms (molecules) while cooling down. Shull and van Steenberg (1985) gave fits 
for the fraction of photo electrons above lOOeV which go into heat, ionization, and excitation respectively. We 
however find that the total energy of all photo electrons above lOOeV is only 1.5% of the total energy of all 
the photo electrons produced in photoionization of neutral hydrogen by a typical quasar flux. Hence we do not 
believe those photo electrons to be significant. The ones in the range of hv t h to lOOeV make up about one third 
of all of the produced photo electrons. Fig. 3 of Shull and Van Steenberg (1985) shows that even for a mostly 
neutral gas, maximally about 30% secondary electrons will go into ionization. We conclude that by leaving out 
secondary ionization we will make an error in the estimates of the ionized fraction of always less than 10%. 

The estimates above were based on the photoionization cross section of hydrogen, but they hold in general 
for the following reasons. First, the cross sections for He and He + scale as z^~ 3 in the high frequency limit, 
and second, nucleosynthesis predicts n# ~ 16n# e . Therefore the photo electrons produced in photoionizations 
of neutral hydrogen are the dominant ones. Clearly, the above discussion only justifies the approximation of 
leaving out the effects of photo electrons, for applications where the radiation spectrum is approximated well 
by a power law as given in equ. 5.1.. 

5.2. Exciting Collisions 

For our intended low density applications the coronal limit (see above) holds and we therefore know that any 
excitation of H, H~, He, He + , and H2 will be followed by a spontaneous decay back into the ground state and 
so ensure that practically every atom or ion will be in its ground state most of the time. An excitation does not 
change the abundance of our model species and does not enter the reaction network explicitly. We, however, 
will consider the dominant ones in the energy balance through appropriate cooling functions. 

5.3. Collisions of Electrons with Atoms and Ions 

Electrons can photo-attach with, ionize, and excite, neutral hydrogen and helium and also ionize, recombine 
with, and excite their ions. We included all radiative recombination and collisional ionization reactions (procs. 1- 
6, 14) except the double ionization of helium since it does not significantly influence the ionization balance 
due to its high threshold. The dominant excitations (see above) are treated through the appropriate cooling 
functions. We also included the photo-attachment reaction of neutral hydrogen (proc. 7). No reference for 
photo-attachment of neutral helium could be found but one can be confident that its effects are by far dominated 
by H~ which is the main species for the formation of molecular hydrogen and thereby controlling the cooling 
of primordial gas at low temperatures. 
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5.4. Proton Collisions with Neutral Helium and Hydrogen 

Below temperatures of ~ leV the neutral atoms are the most abundant species, whereas free protons and 
He + are, even in the non-equilibrium case, more than two orders of magnitude less abundant than H, and He 
respectively (Shapiro and Kang 1987). For the equilibrium case this difference is far more drastic. 

We checked all proton collisions listed in Janev et al. 1987 and find that the rate coefficients typically 
drop drastically for temperatures below lOeV. This is especially true for collisional ionization by protons. The 
exception to the rule are charge exchange reactions as, 

H(n) +p -> p + H(n), n = 1,2,3, .... 

Since we expect nearly all hydrogen atoms to be in the ground state, the n = 1 case will be the most probable 
process. Although it does not enter the reaction network, since it "produces what it destroys", it is an important 
reaction theoretically since it ensures tight thermal and spatial coupling between the protons and the ions due 
to its relatively high rate coefficient of ~ 10~ 8 cm 3 s -1 . 

All rate coefficients for the excitation of neutral hydrogen and helium by protons out of the ground state 
are smaller by many orders of magnitude than the one for excitation by electrons. Since charge neutrality 
requires that there are about as many free electrons as protons, we know that excitation by electrons dominates 
excitations by protons. 

5.5. Collisions of Hydrogen with Helium and their Ions 

An interesting feature of the presented reaction network is that all species which are build up by hydrogen 
nuclei couple with the species formed by helium nuclei only by the free electron fraction. This is due to the 
following reasons; 

1. Electrons are much more effective in ionizing than any other species since their mean velocity is 
(im^/m,;) 1 / 2 ~ 43A 1 / 2 times higher than for an ion with A times the proton mass. 

2. Charge Exchange between H and He + has a rate coefficient ~ 1.9 x 10~ 15 cm 3 s _1 and is according to 
Couchman (1985) negligible. 

3. H~, H^, and H2 are not capable of ionizing helium. They are relatively fragile and have very small 
abundances at temperatures above ~leV where they have kinetic energies far below ionization thresholds. 

A discussion why the reaction network does not include processes where helium or its ions destroy H~, H^~, or 
H2 can be found in the subsequent sections. 

5.6. Neglected Dissociating Reactions of H 

Electron detachment of H~ through neutral helium has according to the ALADDIN database of the IAEADS, 
a cross section of the order ~ 10 -16 cm 2 and is therefore by many orders smaller than the cross section 
(~ 1(T 9 - 1(T 7 cm 2 ; Janev et al. 1987, 7.3.2a) of reaction 

H~ + H -» H 2 *(Z g ) -» H + H + e~. (13) 

Exactly the same argument holds for electron detachment of H~ through molecular hydrogen which has a cross 
section ~ 1(T 18 - l(T 16 cm 2 (ALADDIN database). 



13 



We did not include mutual neutralization of H with He + . The available data from the IAEADS via the 
ALADDIN database, which is only valid for energies above 42eV, indicates that the cross section of mutual 
neutralization of with He + is higher than for mutual neutralization with protons. That is what one 
might naively expect when considering that the Coulomb attraction, acts as the driving force for this reaction, 
competing with the kinetic motion of the reactants. Since He + is heavier than H its mean velocity will be slower 
but the Coulomb attraction is the same therefore mutual neutralization with He + should occur with a higher 
rate. We only leave out this reaction due to the lack of reliable data. However, the expected error is small since 
at low temperatures most of the H~ will react with neutral hydrogen and form molecular hydrogen (proc. 7) 
and at higher temperatures (>8,000K) collisional detachment by free electrons (proc. 14) will dominate the 
destruction of H . 

Neutralization between H~ and He ++ can be neglected due to the high temperature threshhold for He ++ 
formation. Whenever He ++ is formed the temperature is so high that nearly all H~ will be destroyed. 

5.7. Omitted reactions for H2 formation 

Molecular hydrogen formation by excited atom radiative association 

H(n = 2) + H(n = 1) -» H 2 + hu 

has a very small rate coefficient, due to the zero dipole moment of molecular hydrogen (Latter and Black 
1991). Also, for the cosmological applications we are interested in, the n = 2 population is extremely small in 
all cosmological applications making this process less significant than the dominant H2 formation mechanisms 
(procs. 8,10). 

We also neglect the formation path where Hg" first gets formed by excited atom radiative association and 
then forms molecular hydrogen through the charge exchange with neutral hydrogen (proc. 10) (Rawlings et al. 
1993). 

H(n = 2) + H(n = l) -> + e", 

H+ + H(n = l) H 2 + H+. 

This process can in certain circumstances, e.g. dense and hot gas in circumstellar environment, dominate 
formation by excited atom radiative association In our intended cosmological applications these processes, 
however, will never be significant, since the population of excited states is always negligible. 

5.8. Collisional Dissociation of H 2 and 

Since the hydrogen molecule is the main coolant for gas at temperatures below ~leV we pay special attention 
to its destruction mechanisms. We follow here Shapiro and Kang (1987) and assume that at every dissociation 
4.476eV are lost. In the following we will discuss dissociating reactions which we find to be negligible compared 
to the three dissociating reactions we have included in our model and which are illustrated in Figure 2. 

5.8.1. Dissociation of molecular Hydrogen by e~ 
i) + e~ -» (HJ)* -> H" + H 

Wadehra and Bardsley (1978) showed that this dissociative attachment reaction depends strongly on the 
vibrational and rotational states. For low density gas (njj < 10 5 cm~ 3 ) in which nearly all hydrogen molecules 
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will be in their ground state its rate is of the order 10 _15 cm 3 s™ 1 for an electron temperature of leV. For electron 
temperatures below that this rate drops drastically even further. After comparing this to Figure 2 we conclude 
that this process will never play an important role for the destruction of H2 in our applications. It is, however, a 
crucial process in situations with either densities above 10 4 cm -3 or intense ultraviolet radiation fluxes, because 
under those circumstances the excited states will become populated (for vibrational excitation of molecular 
hydrogen in intense UV fluxes see Shull 1978). 

ii) H 2 + e" -> H* + H* + e" , 

Our process 12 is, 

H 2 + e" ->H(ls) + H(ls)+e-, (14) 

and dominates all other reactions of the type, 

H 2 + e- ^H* + H* + e". (15) 

Intuitively one would have expected this because in proc. 15 the incident energy has to be enough to excite the 
produced H atoms. However, the 2s level already lies 10.2eV above the ground state which only few electrons 
have at the low temperature (< leV) where molecular hydrogen exists. 
From the graphs in Janev et al. 1987, (2.2.5) it is evident that also 

/ b'K \ 

H 2 + e~ -► tf 2 - a 3 S+ -» H(ls) + H(ls) + e~ 
V c 3 U u J 

will be dominated by the above reaction 14. Please note that the rate coefficient of reaction 2.2.6 in Janev et 
al. 1987 is a factor 10 to high! (Personal communication with Bill Langer, 1995) 



5. 8. 2. Dissociation of Molecular Hydrogen by H + 

The positive hydrogen ion is able to destroy molecular hydrogen through the charge exchange reaction (11). 
The direct collisional dissociation, 

H+ + H 2 -» 2H +H+, 

has often been left out by former studies. We also were not able to find a rate coefficient for this reaction. 
However, one does not believe this to be important since 

1. at low temperatures where a significant H 2 fraction exists there will be nearly no free protons due to the 
high recombination rate at low temperatures. 

2. dissociation by electrons is likely to be more effective due to their higher velocities in the single temperature 
fluid. 



5. 8. 3. Dissociation of Molecular Hydrogen by He 

Dove and Mandy (1986b) found that He in comparison to H is very inefficient in dissociating H 2 (0, 0). They 
further state that, evidently the collision between two closed shell species tends to cause the molecular bond to 
stiffen, whereas a collision with an atom weakens and loosens the molecular bond. In a follow up paper Dove et 
al. (1987) applied their result to interstellar densities and found that the dissociation through neutral helium is 
negligible and that H 2 -He collisions do not excite the vib./rot. populations of H 2 in low density gas. Following 
thess results we do not include dissociation of H 2 (0,0) by He. 
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5.8.4- Dissociation of Molecular Hydrogen by He + 



The number fraction of He + in primordial gas is about 10%. Due to the high ionization threshold of He + 
of 54.4eV it will not be abundant at low temperatures. Furthermore, from the scaling relations of the rate 
coefficients given in Janev et al. (1987), we know the dissociation of H2 by He + occurs 0.35 times less than the 
corresponding proton reaction. Therefore we can safely conclude that this reaction can be omitted. 

5. 8. 5. Dissociation of Molecular Hydrogen by Hi 

The dissociation rate coefficient for H2 is of the same order and temperature dependence as the one for 
dissociation by neutral hydrogen (Lepp and Shull, 1983). However we do not include this reaction since the H2 
to H fraction (and hence the rate) is always smaller than ~ 10~ 2 and therefore negligible. 

The study of Dove et al. (1987a, b) showed that, in the low-density limit, dissociation of molecular hydrogen 
by collisions with helium at temperatures in the range 2000 - 10 4 K is negligible compared to the H2 - H 
collisions. Only for temperatures close to 10 4 K do the two rates become of the same order. We can, however, 
still neglect this process since helium has a much smaller number density than hydrogen, and second we also 
expect dissociation via the charge exchange reaction with H + to dominate at such high temperatures. 

We do not know the rate for direct collisional dissociation of molecular hydrogen by protons. However, by 
simply looking at the rate coefficients for ionization of H2 and dissociation of H^ by protons (Janev et al.1987, 
3.2.5, 3.2.6) one sees that those rate coefficients drop drastically for temperatures below lOeV. Therefore, we are 
confident that dissociation by protons is negligible. For dissociation by He ++ ions exactly the same arguments 
used for protons are valid (see Janev et al. 1987, 5.2.1, 5.2.3). In addition the charge exchange reaction between 
H2 and He ++ is highly inefficient for temperatures below lOOeV. 

5.9. Collisional induced Absorption 

Molecular hydrogen does not have an electric dipole moment in its ground electronic state. As a result ab- 
sorption of photons can only take place via electric quadrupole transitions and low-density molecular hydrogen 
gas is essentially transparant throughout the visible and infrared part of the spectrum. However, each time 
a collision between two particles occurs, the interacting pair (H2 - H2 , H2 - He, H2 - H) acts as a "virtual 
super molecule" which due to its nonzero electric dipole, can absorb photons with a probability which is much 
higher than that of an isolated H2 molecule. However, this process can, according to the results of Lenzuni, 
Chernoff, and Salpeter (1991), only contribute significantly to the opacity for number densities greater than 
~ 10 18 cm -3 , and hence is negligible at the low densities considered here. 

6. The Cooling and Heating Functions for Optical Thin Gas 

To model the cooling behaviour of the gas correctly, we incorporate the following cooling and heating mech- 
anisms, 

• Compton cooling 

• Recombination cooling due to hydrogen and helium 

• Line cooling of hydrogen and helium 
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• Bremsstrahlung cooling 

• Molecular formation and line cooling 

• Photoionization heating 

• Photodissociation heating 

The cooling functions for Compton, recombination, line, and bremsstrahlung cooling we have used are given in 
AZAN96. The molecular cooling and heating rates are somewhat controversial in the literature and, therefore, 
deserve special attention. 

6.1. Molecular Cooling 
Formation cooling of molecular Hydrogen 

Most of the reaction enthalpy (1.83 eV for process 10 and 3.53 eV for process 8) goes rather into ro- vibrational 
excitation of the molecule than to the kinetic energy of the outgoing particles. We assume that all the reaction 
energy goes into excitation and gets radiated away through subsequent spontaneous decay to the ground state 
(Shapiro and Kang 1987), with the following cooling formula: 

A# 2 Formation = n H (3.53n H -k 8 + 1.38ra H +) eV cuT 3 s _1 

These reactions, however, are important heat sources at high densities (n > 10 8 1/cm 3 ) where collisional de- 
excitation can transform most of the excitation energy into kinetic (thermal) energy. The density dependent 
heating rates are given by Hollenbach and McKee (1979, hereafter HM). 

Line Cooling of Molecular Hydrogen 

Recently Martin, Schwarz, &: Mandy (1996, MSM96 herafter) have derived the molecular hydrogen cooling 
function for H2 - H collisions with a complete set of rate coefficients. Their result differs substantially from 
former cooling functions (e.g. Hollenbach and McKee 1979; Lepp and Shull 1983). Especially at number 
densities eceeding 10 cm~ 3 Lepp and Shull's expression seems to overestimates (underestimate) the cooling at 
low (high) temperatures by an order of magnitude. 

Although, the rate coefficients of MSM96 are only accurate for temperatures above 600 K, some of them have 
been found to agree with QCT calculations at 300 K. For lower temperatures, however, the MSM96 cooling func- 
tion has to be understood as lower limit (personal communication with Peter Martin, 1996). A FORTRAN77 
routine that computes their fitting formula (http : //zeus . ncsa . uiuc . edu : 8080/ abel/PGas/cool . html) 
can be found on our WWW-Site. From Appendix B. we know that the number density of ortho-H2 n£ can be 
very small for applications where our model is valid. Therefore, if one were to estimate the cooling behavior to 
a high accuracy, a separate treatment of the ortho and para-H2 states would be required. 

Photodissociation Heating 

Our model includes the photodissociation processes of H2 , processes (27) and (28). Black and Dalgarno 

(1977) argue that, for typical radiation fluxes, the photodissociation by the two-step Solomon process (27) 
yields 0.4 eV per atom pair. Hence the corresponding heating function is 

T 27 = 6.4 x 10~ 13 n H2 fc 2 7 erg cm" 3 s _1 . (16) 

For direct dissociation into the continua of the Lyman and Werner systems, the situation is similar to 
photoionization heating. Since the reaction channel leading to the excitation of the produced hydrogen atoms 
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is only accessible for excess photon energies greater than 10.2 eV, we can safely assume that all the excess 
energy of the dissociating photon will be shared as kinetic energies of the produced hydrogen atom pair. Hence 
the heating functions can be written, 



where i(v) denotes the radiation intensity and the integral is evaluated for the ortho and para-H2 component 
separately because of their different threshold energies. 



In this section we will use data from a Id high resolution study of a cosmological sheet with wavelength 
4Mpc, to discuss non-equilibrium effects and the chemical dynamics in primordial gas. The data is taken from 
Anninos & Norman 1996. 

Very massive pancakes have strong virialization shocks (y s ~ lOOkm/s), which leads to the destruction of 
primordial pre-shock hydrogen molecules. The post shock gas cools faster than it recombines and leaves the gas 
with a significant free electron fraction even at low temperatures. This can clearly be seen in Figure 3 where 
the fractional abundances of all species in the post shock gas are plotted as a function of temperature. There 
we also see that is always much less abundant than H~. The doubly negative helium ion recombines very 
fast and therefore shows the non-equilibrium effect at somewhat lower temperatures than hydrogen. These 
abundances agree well with the ones given in Shapiro and Kang (1987). 



The right hand side of an individual rate equation is, for collisional processes, given by a sum of terms kijiiiiij, 
where kij is the rate coefficient for the reaction between species i and j, which is taken positive if it produces 
the species under consideration, or negative if it destroys it. To see what processes dominate the chemistry of 
species I, let us consider the rate kijUiUj divided by ni, which we call rate per /-atom (molecule). Comparing 
this to the sum over all I producing/destroying processes (the evaluated right hand side) much insight about 
the ongoing processes can be gained. Actually this quantity equals the inverse reaction time scale and hence 
provides a measure of the time in which a particular species gets into equilibrium. We show such plots for 
H2 , , and H~ (Fig. 4, 6, and 5). The complete H2 chemistry in our model is illustrated in Fig. 4. All 
molecular hydrogen producing reactions (8, 10, and 19) are shown as solid lines as well as their sum (labeled 
formation). The destroying processes (11, 12, and 13) are depicted by dashed - dotted lines. It is clear that the 
molecular hydrogen abundance at temperatures above ~ 7000 K is a result of the balance between producing 
and destroying processes, what one could call a "collisional equilibrium". It can also be seen that from 0.5eV to 
leV the charge exchange reaction of H2 and free protons is the most efficient molecular hydrogen destruction 
path. At higher temperatures, molecular hydrogen are destroyed most efficiently by free electrons. Obviously 
here in the pancake collapse, destruction by neutral hydrogen atoms has a negligible influence on the molecular 
hydrogen abundance, which is certainly due to the high (non-equilibrium) abundance of free protons which will 
not be found as drastically enhanced, compared to the equilibrium abundances, for weak shock waves. For the 
production of H2 we find that proc. 8 and 19, are always dominated by the fast H~ formation path (proc. 8). 

In summary, we find that the molecular hydrogen fraction is governed dominantly by three processes. Forma- 
tion via H , destruction through charge exchange with free protons, and destruction by free electrons. These 
insights can be used also for analytical estimates of the molecular hydrogen fraction during pancake collapse. 
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7. Application and Discussion of the Model 



7.1. H 2 Chemistry 
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Let us illustrate this briefly on the molecular hydrogen abundance for temperatures, 7000 K < T < 10 4 K. 
Here the H2 abundance can obviously be derived through the equilibrium condition dnn 2 /dt = with, 

—— ~ k 8 n H -n H - k n n H2 n H + = 0, =>- n H2 - ~, 

dt fen x 

where x denotes the ionized fraction. Looking up the rate coefficients we can read off the qualitative tem- 
perature dependence and find that with a decreasing temperature in the range 6000 K < T < 2 x 10 4 K, the 
molecular hydrogen abundance increases. For lower temperatures, collisional processes are not efficient enough 
in destroying molecular hydrogen and it will get produced as long as there exists a significant amount of H~ 
ions at the level of a few percent. 



7.2. The H and Chemistry 



The H~ ion is found to get into chemical equilibrium faster than the dynamical time scales (Courant, free- 
fall, and cosmological times). Hence the overall destruction and production lines in Fig. 5 are right on top 
of each other. Our model includes one H producing reaction, the radiative attachment (proc. 7), and six 
destroying reactions. Fig. 5 clearly illustrates that the processes 17 and 19 play no role, processes 15 and 16, 
little role, and processes 7, 8, and 14 the main role, for the H - chemistry. Since the atomic data we used 
for the main processes 7, 8, and 14 are very accurate at all temperatures we ensured that the model will give 
also accurate predictions for the H~ abundance. An analogous plot for the H^ chemistry is given with Fig. 
6. Obviously H^~ has been in chemical equilibrium for all temperatures in our pancake study. All of its three 
production mechanisms contribute in individual temperature regimes. The radiative association reaction proc. 
9 contributes strongly for temperatures below a few thousand Kelvin and dominates at higher temperatures as 
~ 2 x 10 4 K. For the remaining interval charge exchange between H2 and H + is the dominant H^ producing 
reaction. Furthermore it is clear from Fig. 6 that one can safely neglect the destroying mechanism proc. 
19, because it does not contribute significantly, at any temperature, to the overall destruction of H^ . It is 
also evident that at all temperatures below ~ 8 x 10 3 K nearly every H^ molecule that gets produced will 
be converted to a neutral hydrogen molecule by proc. 10. At higher temperatures H^ will be destroyed by 
electrons. The data of Schneider et al. for the rate coefficient of this process, and consequently our fit, is quite 
inaccurate at such high temperatures and we do not expect the H^ abundance at temperatures higher than 
10 4 K to be reliable. We see that process 19 was never important in determining the abundances of H~, H^ , 
and of H2, and can therefore be neglected in future applications. 



7.3. A Minimal Model 



From the preceding section we clearly see that not all reactions are equally important. We devised a minimal 
model that incorporates only the essential processes need to accurately model the formation of molecular 
hydrogen. From Figure 4 it is evident that the only reaction involving H^ which contributes strongly to the 
H2 chemistry is the charge exchange reaction (11). Its rate is independent of the H^~ abundance since it is a 
H2" producing reaction and we can conclude that at least for strong shocks the H2 abundance is independent 
of the H^ chemistry. For a minimal model we can use this and leave out reactions 9, 10, 17, and 19 so as to 
avoid solving a H^ rate equation. Further by looking at Figure 5 we find the process (15) to be negligible and 
process (16) only to be marginally important. Clearly (13) can be left out as well since it obviously does not 
contribute to the H2 chemistry. Furthermore, we note that it is not necessary to include the contributions from 
the reaction involving H and H2 to the major species H, He, and their ions. With these simplifications we 
eliminated 7 of the 19 reactions and reduced the reaction network to the following, 
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dn 



H 



dt 
dn H + 

dt 
dn Hc 

dt 
dn Re + 
dt 
dn He ++ 
dt 

dnn 2 



/c 2 n H +nc - fcinH«e (18) 

kinnn e - /c 2 n H +n e (19) 

k4n Hc +n e - kznuene (20) 

k3n-n e n e + ken Re ++n c - kAn Ke +n e (21) 

k$n ne +n e - k6n He ++n e (22) 

ksn H -nn-nn 2 (kiin u + + ki 2 n e ) , (23) 



dt 

where the number density of H~ is given by the equilibrium condition 

n _ = k 7 nun c 

fc 8 n H + fei 6 n H + + kun e ' 

To derive the number denisty of H~ by its equilibrium value is justified by comparing the figures, 4 and 5 . 
Since they show the inverse of the reaction time scales of the different processes it is clearly evident that all 
reactions determining the H~ abundance occur on much faster time scales than those responsiblefor the H2 
chemistry. We have checked this minimal model extensively and find the H2 abundance to generally aggree 
with the full model to within a few percent over the entire temperature of the individaul application. 



8. Conclusions 

We have derived a reliable time dependent chemistry and cooling model for primordial gas that, in connection 
with our new numerical method (AZAN96), proves as a powerful tool to investigate primordial structure 
formation in multidimensional numerical calculations. The model is designed to be applicable for densities 
below 10 4 cm -3 and temperatures < 10 8 K. We have discussed the influence of the orth-H2 to para-H2 ratio on 
the cooling behaviour of the gas, derived new fits to molecular data, and developed a minimal model capable 
of describing primordial gas in applications where no external radiation fields are considered. 



9. Related Websites 

At our WWW-Site, The LC A Cooling Model (http : //zeus . ncsa . uiuc . edu : 8080/ abel/PGas/LCA-CM . html), 

we present all the used atomic and molecular data in great detail. Further we provide there FORTRAN routines 
that compute the rate coefficients and sove the rate equations. Tom Abel continues to collect molecular data re- 
lated to primordial gas at The Primordial Gas Chemistry Page (http: //zeus .ncsa. uciu. edu: 8080/ abel/PGas/). 
Further Dima Verners Atomic Data Page (http://www.pa.uky.edu/ verner/atom.html) represents a 
superb reference for atomic data. 

We would like to thank Dimitri Mihalas, Mordecai-Mark MacLow, Zoltan Haiman, Max Tegmark, and Eve- 
lyne Roueff for helpful discussion. Furthermore, we are grateful for useful correspondence with Bill Langer, Phil 
Stancil, Jonathan Rawlings, Alex Dalgarno, Peter Martin, and Stephen Lepp. This work was done partly under 
the auspices of the Grand Challenge Cosmology Consortium funded by NSF grant ASC-9318185. Tom Abel 
happily acknowledges the hospitality of the Max Planck Institut fur Astrophysik and the Insituto Astrofisica 
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A. Reactions and Rates 



Here we present all the rate coefficients used in our model. All fits are accurate for temperatures ranging 
from IK to 10 s K. Since we are interested in numerical applications we did pay more attention to the accuracy 
of the fits than to the simplicity of the formulas. A FORTRAN program that computes these rate 
coefficients (http : //zeus .ncsa.uiuc . edu: 8080/ abel/PGas/LCA-CM.html) can be obtained at our WWW 
site. The temperatures are in eV unless stated otherwise. 



Table 1: Collisional Ionization and Radiative Recombination of Hydrogen and Helium. 

(1) H + e~ -> H+ + 2e~ Janev et al. 1987 (2.1.5) 

fcj = cxp[-32. 71396786 + 13.536556 ln(T) - 5.73932875 ln(T) 2 + 1.56315498 ln(T) 3 - 
0.2877056 ln(T) 4 + 3.48255977 x 10~ 2 x ln(T) 5 - 2.63197617 x 10~ 3 x ln(T) fi + 
1.11954395 x 10~ 4 ln(T) 7 - 2.03914985 x 10 -6 ln(T) 8 ] cm 3 s _1 . 

(2) H+ + c~ -> H + 7 Our fit to data given by Ferland et al. (1992) 

k 2 = cxp[-28. 6130338 - 0.72411256 ln(T) - 2.02604473 x 10~ 2 ln(T) 2 - 

2.38086188 x 10~ 3 ln(T) 3 - 3.21260521 x 10~ 4 ln(T) 4 - 1.42150291 x 10~ 5 ln(T) 5 + 
4.98910892 x 10~ 6 ln(T) 6 + 5.75561414 x 10~ 7 ln(T) 7 - 1.85676704 x 10~ 8 ln(T) 8 - 
3.07113524 X 10" 9 ln(T) 9 ] cm 3 s _1 . 

(3) He + c~ -> Hc+ + 2e~ Janev et al. 1987 (2.3.9) 

k 3 = cxp[(-44. 09864886 + 23.91596563 ln(T) - 10.7532302 ln(T) 2 + 3.05803875 ln(T) 3 - 
0.56851189 ln(T) 4 + 6.79539123 x 10~ 2 ln(T) 5 - 5.00905610 x 10~ 3 ln(T) 6 + 
2.06723616 x 10~ 4 ln(T) 7 - 3.64916141 x 10~ 6 ln(T) 8 ) cm 3 s _1 . 

(4) Hc+ + c~ He + 7 Cen (1992) and Aldrovandi & Pequignot (1973) 

Radiative : k 4r = 3.925 x i - 13 T-°' 6353 cm 3 s" 1 

~ , , q m 3 ( 48.596 cV\ r /8.10eV\ 
Dielectronic : k 4d = 1.544 x 10~ 9 T" T exp ( J 0.3 + cxp ( J 

(5) Hc+ + c~ -> Hc++ + 2e~ Aladdin Database (1989) 

k 5 = cxp[-68. 71040990 + 43.93347633 ln(T) - 18.4806699 ln(T) 2 + 4.70162649 ln(T) 3 - 
0.76924663 ln(T) 4 + 8.113042 x 10~ 2 ln(T) 5 - 5.32402063 x 10~ 3 ln(T) 6 + 
1.97570531 x 10~ 4 ln(T) 7 - 3.16558106 x 10~ 6 ln(T) 8 ] cm 3 s _1 . 

(6) Hc++ + e~ -> He+ + 7 Cen (1992) 

1 / T \-0.2 / / T \°' 7 \ _1 

fe 6 = 3.36 x 10" 10 T-2 ( 1 + _ cm 3 s -1 , where T is in K. 

VlOOOK/ ^ Vl0 f 'K/ J 
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Table 2: The formation paths of H2. 

(7) H + c- -> H~ + 7 



This work from data by Wishart (1979) 



k 7 



1.429 X io- 1 8t - 7620 T - 15231o «io( t )T- 3 - 274x10 2 lo Sio( T ) cm 3 s^ 1 for T < 6000 K 

3.802 X 10- 17 T°- 19981 °Sio( T )dex (4.0415 X 10~ 5 logf (T) - 5.447 X 10~ 3 log 4 (T)) cm 3 s" 1 otherwise. 



(8) H + H — » H2 + e Our Integration of data from Janev et al. (1987) 

T > 0.1 cV : fc 8 = exp[-20.06913897 + 0.22898 ln(T) + 3.5998377 x 10~ 2 ln(T) 2 - 
4.55512 x 10~ 3 ln(T) 3 - 3.10511544 x 10~ 4 ln(T) 4 + 1.0732940 x 10~ 4 ln(T) 5 - 
8.36671960 x 10~ 6 ln(T) 6 + 2.23830623 x 10~ 7 hi(T) 7 ] cm 3 s _1 . 

T < 0.1 cV : fc 8 = 1.428 X 10~ 9 cm 3 s" 1 . 

(9) H + H+ ->■ H+ + 7 Shapiro & Rang (1987) 

f 3.833 x 10~ 16 x T 18 cm 3 s _1 , for T < 0.577 eV, 

9 ~ \ 5.81 x 10" 16 x (0.20651 T)-°- 2891xl °s(o.2065ixT) cm 3 s -i 3 for T > 577 eV . 

(10) + H -> H 2 + H+ Karpas ct al. (1979) 

fclO = (6.4 ± 1.2) X 10~ 10 cm 3 s _1 . 
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Table 3: Other Collisional Processes involving H , , and H2. 

(11) H 2 + H+ -> H+ + H This work 

ln(fcn) = -24.24914687 + 3.40082444 ln(T) - 3.89800396 ln(T) 2 + 2.04558782 ln(T) 3 - 
0.541618285 ln(T) 4 + 8.41077503 x 10~ 2 ln(T) B - 7.87902615 x 10~ 3 ln(T) 6 + 
4.13839842 x 10~ 4 m(T) 7 - 9.36345888 x 10~ 6 ln(T) 8 cm 3 s _1 . 



(12) H 2 + e~ -» 2H + e~ Donahue and Shull (1991) 

n i , 102.124 K, , , 
fci2 = 5.6 X 10 _11 T2 exp( - ) cm 3 s i , where T is in K. 



(13) H 2 + H -> 3H Dove and Mandy (1986) 

fci3 = 1.067 x io- 10 T 2 012 x exp(-(4.463/T)(l +0.2472 T) 3 ' 512 ) cm 3 s" 1 . 

(14) H~ + c~ -> H + 2 er Jancv et al. (1987, 7.1.1) 

k 14 = cxp[-18.01849334 + 2.3608522 ln(T) - 0.28274430 ln(T) 2 + 1.62331664 x 10~ 2 ln(T) 3 - 
3.36501203 x 10~ 2 ln(T) 4 + 1.17832978 x 10~ 2 ln(T) 5 - 1.65619470 x 10~ 3 ln(T) 6 + 
1.06827520 x 10~ 4 ln(T) 7 - 2.63128581 x 10~ 6 ln(T) 8 ) cm 3 s -1 . 

(15) H~ + H — > 2H + c~ Our modification to the Janev et al. (1987) data 

T > 0.1 cV : feiB = exp[-20.37260896 + 1.13944933 ln(T) - 0.14210135 ln(T) 2 + 

8.4644554 x 10~ 3 ln(T) 3 - 1.4327641 x 10~ 3 ln(T) 4 + 2.0122503 x 10~ 4 ln(T) 5 + 
8.6639632 x 10" 5 ln(T) 6 - 2.5850097 x 10~ s hi(T) 7 + 2.4555012 x 10~ 6 ln(T) s - 
8.0683825 x 10~ s ln(T) 9 ] cm 3 s" 1 , 

T<0.1cV : fci 5 = 2.5634 X 10" 9 X T 1 ' 78186 cm 3 s _1 . 

(16) H~ + H+ -> 2H Dalgarno and Lepp (1987) 

k 16 = 7 x 10~ 8 ( ] 2 cm 3 s~\ where T is in K. 

VlOOK/ 

(17) H" + H+ -> H+ + e" Our modification to fit of Shapiro and Kang (1987) 



2.291 x io- 10 T-°- 4 cm 3 s" 1 , for T < 1.719 eV 

8.4258 X 10- 10 T- 1 - 4 exp(-1.301/T) cm 3 s" 1 , otherwise. 



(18) H+ + e~ -» 2H Our fit to the Schneider et al. (1994) data 

_ f 1.0 x 10~ 8 cm 3 s" 1 , for T < 617 K, 

18 ~ \ 1.32 X 10- 6 T- (L76 cm 3 s" 1 , for T > 617 K, where T is in K. 

(19) H+ + H" -> H 2 + H Dalgarno and Lepp (1987) 



fcig = 5 x 10 



^ip^Ky 
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Table 4: Photoionization and Photodissociation. 

(20) H + 7 -> H+ + e~ , 

(22) Hc+ + 7 -* He++ + e~ Ostcrbrock (1974) 

A) ( v \ 4 cxp[4 — 4(arctane)/e] ^ ,„_is 2 

<720,22 = — T -, n , [ J , where A = 6.30 X 10 18 cm 2 , 

Z 2 \VthJ 1 - cxp(-27r/e) 

e = V^Mft ~ !> hv th = 13.6Z 2 eV. 

(21) He + 7 -> Hc+ + e~ Osterbrock 1974 



/ / x -2.05 / x -3.05\ 
0-21 = 7.42 X 10~ 18 I 1.66 ^ J - 0.66 ^ J j cm 2 , for 1/ > v th . 

(23) H~ + 7 -> H + e~ Shapiro and Kang (1987) 



3 1 

(T23 = 7.928 x 10 5 (^ - v th )? — cm 2 , for hv > /u/ t/l = 0.755 eV. 
v 6 



(24) H 2 + 7 -> H+ + e- O'Neil and Reinhardt (1978) 

{0, for hv < 15.42 eV, 

6.2 X lQ- ls hv - 9.4 X 10~ 17 cm 2 , for 15.42 eV < hv < 16.50 eV, 

1.4 x 10- is hu - 1.48 x 10~ 17 cm 2 , for 16.50 eV < hv < 17.7 eV, 

2.5 X 10- 14 (/n/T 2 ' 71 cm 2 , for /if > 17.7 eV. 

(25) H+ + 7 -> H + H+ Our fit to Stancil (1994) 

l°gio( CT 25 S ) = (-1-6547717 x 10 6 + 1.8660333 x 10 5 ln(i/) - 7.8986431 x 10 3 ln(^) 2 + 
148.73693 ln(i/) 3 - 1.0513032 ln(z/) 4 ) 

(26) H+ + 7 -> 2H+ + e - Shapiro and Kang (1987) 

log 10 (o- 2 6) = -16.926 - 4.528 x 10~ 2 hv + 2.238 x 10~ 4 (/«/) 2 + 4.245 x 10~ r (hv) :i , 
for 30 eV < hv < 90 eV. 

(27) H 2 + 7 ^ H* ^ H + H This work 

Neglecting sclfshiclding, fc27 = 1.1 X 10 8 j(F) s _1 , 

where j(v), is the radiationflux in ergs s _1 cm -2 at hv = 12.87 eV. 

(28) H 2 + 7 -> H + H This work 

CT 28 - — (^28 + °28 ) + l 1 — r)( CT 28 + CT 28 ) 

2/ + 1 2/ + 1 
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B. Ortho-H2 to para-H2 Ratio 



Many of the collisional rates have only been computed for the para configuration and are applicable only 
if para-H 2 is is the dominant modification (Flower and Watt 1984, hereafter FW). For our purposes it plays 
another significant role in that its selection criteria for allowed rotational transitions are different. This is crucial 
to determine the exact cooling behavior at low temperatures, since low density primordial gas will cool mostly 
in the rotational lines E% — E\ = 844. 5K for the ortho configuration and E2 - Eq = 509. 9K for para-H2 . At 
formation the two different modifications will be abundant corresponding to their statistical weights, namely 

njortho) = n(J = 1) = g 
n(para) n(J = 0) 

FW argue for a single interconversion mechanism at low densities, 

H 2 (ortho) + H+ -» H 2 (para) + H+ + 170.5 K, (B2) 

and suggest a temperature independent rate coefficient of, k op = 3 x 10~ 10 cm 3 s _1 . If the J = 2 state would 
be significantly populated, para-H 2 could be transformed to ortho by, 

H 2 (J = 2) + H+ -» H 2 (J = 1) + H+ + 170.5 K. 

However, at the low density limit mostly the ground level will be populated and equ. (B2) will determine the 
ortho to para-H 2 ratio. 

Assuming that we have a constant total number density of molecular hydrogen, n 2 , the rate equation for 
para-H 2 , n 2 simply becomes, 



^ = n 2 nxfc op , (B3) 



where n 2 , n, x, denote the ortho-H 2 number density, the neutral hydrogen number density, and the ionized 
fraction respectively. The time evolution of the ionized fraction is determined from 

— = -k 2 nx 2 (B4) 

to be 

x(t) = — 7 . (B5) 

w 1 + x k 2 nt v ; 

Using n 2 = n 2 + n 2 , and equ. (B5), we derive the simple solution for the number density of para-H 2 

n\ = n 2 (l- l x ^o P /k 2x -k ov /k^ =n2 (i_ 3 x Wfc 2 (n Hk2t + 1^) W ^ _ (B6) 

From this we clearly can see that if 

• k op <Sk 2 , ortho-H 2 to para-H 2 will not change within one recombination time. 

• k 2 <^k op , ortho-H 2 to para-H 2 will change dramatically within one recombination time. 
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Since k.2 is well fitted by 1.8 x 10~ 10 T~ 065 cm 3 s _1 , it convincingly shows that always the latter case applies. 
Obviously the recombination time scale also sets the formation time scale of H2 since the electrons act as 
katalysts (see also Abel 1995, and Tegmark et al. 1996). Hence when H2 gets formed by the gas phase reactions 
(7) through (10), it might immediately get converted to its para configuration. The subsequent cooling of the 
gas will therefore happen mostly in the E2 ~ Eq line of para-H2 . For completely neutral hydrogen gas, however, 
at the low density and temperature limit the ortho-H2 to para-H2 ratio will be given by 9exp(— 170.5 K/T) 
(Mandy & Martin 1993). Sun & Dalgarno (1994) computed the rate coefficients for odd transitions of low lying 
rotational levels by collisions with atomic hydrogen, in the temperature range 30 K — 1000 K. For temperatures 
below 300 K these reactions have rate coefficients <C 10~ 17 cm 3 s _1 , i.e. for such low temperatures in low density 
gas (riH ~ 1 cm~ 3 the corresponding reaction time scales are 3> lOGyrs. Hence if the gas can cool faster than 
it recombines below 300 K the above result is unchanged. However, the ortho-H2 to para-H2 will be different in 
different environments. Once the complete set of rate coefficients for collisional excitation and dissociation of 
H2 molecules by H atoms, given by Martin & Mandy (1995), are extended to temperatures below 450 K, and 
extended to H2 - H2 collisions as well as H2 - H + reactions, it will be possible to study the ortho-H2 to para-H2 
ratio in primordial gas for the entire range in which our reaction network is applicable. In summary, we would 
like to stress the need for a fit to the H2 cooling function that depends not only on temperature and density, 
but also on the ortho-H2 to para-H2 ratio. 
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Fig. 1. — The rate coefficients for photo-attachment to H (dashed line), formation (dotted line) and the 
two H2 forming reactions, proc. (8) and (10) (solid lines). 
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temperature in K 



Fig. 2. — Rate coefficients for the dissociation of molecular hydrogen by neutral hydrogen (dot-dashed line), 
and free electrons (dashed line). The light solid line depicts the rate coefficient for charge exchange between 
H + and H2, given by Donahue and Shull (1991), whereas the thick solid line is our numerical integration from 
the Janev et al. data. 
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temperature(K) 



Fig. 3. — The fractional abundance of nine species in a collapsing pancake of wavelength 4 Mpc are shown vs. 
temperature. The data is taken from Anninos and Norman (1996). Clearly the non-equilibrium enhancement 
of free electrons at low temperatures can be seen. Note also the similar non-equilibrium behavior of the He ++ 
fraction. 
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Fig. 4. — The inverse of the reaction time scale for all processes in our model, that involve hydrogen molecules. 
These relative rates for the producing collisional processes, (8), (10), and (19), are illustrated by solid lines. 
The H2 destroying processes, (11), (12), and (13), are shown with dot-dashed lines. The thick solid (thick dot- 
dashed) line depicts the inverse of the sum of reaction time scales of all H2 producing (destroying) processes. 
The divergence of these producing and destroying curves at low temperatures indicates that the H2 molecules 
are out of equilibrium. 
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Fig. 5. — The relative rates are shown for the collisional processes involving H , including the single H~ 
producing reaction (7), and all its destruction processes (8), (14), (15), (16), (17), and (19). It is evident that 
H reached its equilibrium abundance within the hydrodynamical and cooling time scale since the sum of the 
production rates equals the sum of the destruction rates. Here the thick dot-dashed lines illustrate the major 
H destroying processes, (8) and (14). 



34 




Fig. 6. — The inverse of the reaction time scale for the collisional processes involving • The producing 
processes (9), (11), and (17) are depicted by solid lines and the destroying reactions (10, (18), and (19) 
correspond to the dot-dashed lines. Analogous to H~ (figure 5), it is clear that H^" reached its equilibrium 
abundance within the hydro dynamical and cooling time scales. 
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